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Abstract. In this paper we consider three different synchronization problems consisting 
1 in designing a nonlinear feedback unidirectional coupling term for two (possibly chaotic) 

dynamical systems in order to drive the trajectories of one of them, the slave system, 
to a reference trajectory or to a prescribed neighborhood of the reference trajectory of 
the second dynamical system: the master system. If the slave system is chaotic then 
synchronization can be viewed as the control of chaos; namely the coupling term allows 
to suppress the chaotic motion by driving the chaotic system to a prescribed reference 
trajectory. Assuming that the entire vector field representing the velocity of the state 
can be modified, three different methods to define the nonlinear feedback synchronizing 
^ 1 controller are proposed: one for each of the treated problems. These methods are based on 

' results from the small parameter perturbation theory of autonomous systems having a limit 

1 cycle, from nonsmooth analysis and from the singular perturbation theory respectively. 

Simulations to illustrate the effectiveness of the obtained results are also presented. 
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1 Introduction 



^ . , In recent years in the literature on dynamical system analysis a considerable 

' attention has been devoted to the problem of synchronization of coupled non- 

linear dynamical systems (see e.g. [I], [3], [11], [19], [22], [23 )• of the 
most effective methods for solving such problem consists in designing a feed- 
back coupling term which drives the trajectories of one of the two systems 
(the so-called slave system) to a prescribed reference trajectory of the second 
one (named master system). Examples of such approach can be found, for 
instance, in |21j where the coupling term is represented by a linear feedback 
of the tracking error. In [5S] a bidirectional linear coupling term is proposed 
to synchronize two chaotic systems. An approach to synchronization based 
on the classical notion of observers, when the state is not fully available. 
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can be found in [5] and [15] . In many cases, when one deals with nonlinear 
chaotic dynamical systems, the interest is that of steering any trajectory of 
the chaotic system to an equilibrium point or to a limit cycle of the same sys- 
tem or of another coupled system, see [5] , [TD] , [13] and [55] . For an adaptive 
control approach using a linear reference model we refer to [26j . Finally, in 
[12] it is stressed how the carelessness application of the mathematical tools 
of the synchronization theory can lead to incorrect results. 

In this paper, under the condition that the entire vector field of the veloc- 
ity can be modified, we aim at designing a nonlinear feedback unidirectional 
coupling term, based on the state model system, in such a way that all the 
trajectories of the slave system are steered to a prescribed reference trajec- 
tory of the master system. In other words, the coupling term makes stable, in 
a sense that it will be precised for each problem in the following, a prescribed 
trajectory of a dynamical system with respect to the trajectories of an other 
dynamical system by coupling these systems by means of a suitably defined 
nonlinear feedback coupling term. 

Following the linear feedback approach of the previously cited references, 
we provide examples of how it can be possible, by means of different mathe- 
matical theories, to define a nonlinear unidirectional feedback coupling term 
in order to determine a prescribed dynamical behavior to some classes of 
nonlinear dynamical systems in the case when this term can affect each com- 
ponent of the velocity vector field. 

The problems that we will treat in this paper are illustrated in the sequel. 
The first one is the problem of the synchronization of the phase of a limit 
cycle of an autonomous system with that of the limit cycle of the same 
period of another autonomous system. The feedback design is based on 
classical results due to Malkin [14] on the existence of periodic solutions of 
an autonomous system perturbed by a small parameter nonautonomous term 
and on their behavior when the perturbation disappears, namely when the 
parameter tends to zero. Many authors, see for instance [2], [17] and [20] and 
the extensive references therein, have considered the problem of the control of 
the balance between the phases of the subsystems state variables oscillations 
by coupling the subsystems in different ways, i.e. by suitably balancing the 
energy due to the interaction. In particular in [ll a dynamic feedback coupling 
term for phase locking of non identical oscillators is presented. 

The second and third problem consist in the synchronization of the tra- 
jectories to a reference trajectory of the master. To solve these problems 
we adopt two different feedback laws based on a sliding manifold approach. 
First, we define a static discontinuous feedback coupling term with a gain 
depending on the bounded set of the initial conditions for the trajectories 
of the slave system. It is defined by means of the signum function of the 
tracking error, that is by the signum of the difference between a trajectory 
of the slave system and the reference trajectory. By means of a suitably 
defined nondifferentiable Liapunov function and its subdifferentiability prop- 
erties ^ we can prove that any trajectory originating from a given bounded 
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set converges to the reference trajectory in an estimated finite time which 
depends on the feedback gain. It is worth to observe that, since the right 
hand side of the slave system is discontinuous with respect to the state, it 
is necessary to introduce a suitable concept of solution for this system, in 
fact we consider solutions in the sense of Filippov 9 . The discontinuity 
along the reference curve yo = yo{t),t > 0, of the feedback law makes the 
coupled system robust against modelling errors and external disturbances, in 
the sense that the tracking error tends to zero in finite time also in presence 
of modelling imprecision and disturbances, if the gain is sufficiently large. 
Moreover, observe that it is possible to get the reference trajectory at any 
prescribed speed by suitably increasing the gain. Since the implementation of 
the associated switchings across the reference curve is necessarily imperfect, 
in practice switching is not instantaneous and the value of the tracking error 
e{t) = x{t) — yo{t),t > 0, is not perfectly known; this leads to the chattering 
phenomenon which is the main drawback of this feedback law, namely the 
trajectory of the slave system rapidly oscillates around the reference trajec- 
tory. This is a quite undesirable effect, in fact each component of the signum 
of the tracking error in the coupling term switches very fast between the 
positive and negative value of the corresponding gain, which is not a feasible 
behavior for the physical implementation of the control law. In the frame- 
work of synchronization of chaotic system a feedback of this type has been 
used in [28], where the chattering phenomenon has been also emphasized. 

To avoid the undesirable chattering phenomenon we then consider a dy- 
namic feedback, as introduced in [4J, defined by means of a differential equa- 
tion involving the slave system and the reference trajectory. This equation 
depends also on a small parameter and it satisfies, under general conditions, 
all the assumptions of the classical singular perturbation theory on infinite 
intervals (TB]. As we will see this ensures that any trajectory of the slave 
system approaches the reference trajectory within any prescribed error. 

Finally, we present a simulation for each of the considered problems 
which illustrates the effectiveness of the obtained results. Precisely, for the 
first problem we consider, both for the master and slave system, the same 
FitzHugh-Nagumo type equation which has an asymptotically stable limit 
cycle and by implementing our method we synchronize the phase of the slave 
with that of the master. For the second and third problem we have considered 
a chaotic neural network as the slave system 6 , and as master system a neu- 
ral network which possesses a globally asymptotically exponentially stable 
periodic solution which is taken as reference trajectory [18 . 

The paper is organized as follows. In Section [2] we treat the problem 
of the phase synchronization of two self-oscillating nonlinear dynamical sys- 
tems. In Section [3] for two nonautonomous systems we steer any trajectory 
of one of these two systems to a prescribed trajectory of the second one by 
means of a high gain discontinuous feedback of the tracking error. In Section 
E] we consider the same systems and we design a dynamical feedback cou- 
pling term which drives any trajectory of the slave system to any prescribed 
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neighborhood of the reference trajectory. Finally, in Section[5]we present the 
simulations which illustrate the meaning of the obtained results. 



2 Phase synchronization of hmit cycles 

In this section we consider two autonomous systems 

i = f{x), (1) 

where / e C2(M",R"), and 

y = 9iy), (2) 

where g G C^(M", E"). We assume that they have orbitally stable limit cycles 
xo{t) and yo{t), respectively, of the same period T. According to the termi- 
nology adopted in the literature for the problem that we will treat here, we 
will refer to systems (P) and ^ as slave and master system respectively. We 
are interested in the phase synchronization of the slave system to that of the 
master system by adding a coupling term to ([T]). Specifically, we will show 
that for every /i > we can define a function <i>^ : R" x M" — > R" such that 
the system 

fix) + <i>f,{x,yoit)), t>0, 

has a unique asymptotically stable T-periodic solution x^ in a neighborhood 
of {xo{t) : t e [0,T]} satisfying the property 



2 r 2 

kA'('r) - 2/o(t)P dr - min / |.to(t + s) - ?;o(t)|^ rfr 

sela,T] Jo 



< M- (3) 



For this we assume that the Floquet multiplier of the linearized systems 
around XQ{t) and yo{t), equal to 1 is simple and that the others n — 1 are 
inside of the unit open circle. We consider the following system 

x^ f{x) + e{\x~yo{t)f- 
l rT (4) 
se[o,T] T Jq J 

where e and 5 are positive scalar parameters. 
We can prove the following result. 

Theorem 1 Assume, that the equation 



T i-T 



|xo(r + 6l) -yoMrdr = min / |xo(r + s) - j/o(r)|^rfr (5) 

has an unique solution 6q G [0,r]. Then for every fi > there exists (5^ > 
such that for every 6 G [0, (5^] there is > for which the following results 
hold for e G (0,£a). 
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1 ) System ^ possesses a unique asymptotically stable T -periodic solution 
such that 

x^{t) £ Afs^Xxo), for any t e [0,T], 



wtiere J\fs (xn) ~ < x ^ M" : inf \x — xo(t)\ < Sn)- denotes the 5n-neiqhhor- 
" te[o,T] V /I 

hood in R" of the limit cycle xq, 

2) the solution satisfies property 

To prove this theorem we need the following result due to I. G. Malkin [T3], 
which is one of the main tools for the study of the synchronization of coupled 
systems (see [3])- Consider the system 

i = fix) + e-f{t,x) (6) 

where 7 G C"'^(M x M", R"), assume that 7 is T-periodic with respect to time. 
Then it is possible to show that system 

z = -(/'(xo(i)))*^ 

has a T-periodic solution z* such that 

(z*(t),io(t)) = 1, for any t e [0,T]. (7) 

Let us introduce the function F : M ^ M as follows 

F{e)= I {z*{T),-/{T-e,XQ{T)))dT, for any 61 eM. 



We can now formulate the following result. 

Theorem 2 (|14|. Theorems pp. 387 and 392) Assume that for suffi- 
ciently small e > system 0) has a continuous family s x^ of T-periodic 
solutions satisfying the property 

Xe{t) xo{t + Oo) as e^O (8) 

then F{eo) = 0. Moreover, if F{eo) = and F'{0o) 7^ then ^ holds 
true and the solutions x^ are asymptotically stable or unstable according to 
whether F' {60) is negative or positive. 

Proof of Theorem 1. For (5 > let 

/ 1 .T 



js{t,x) = l\x-yo{t)\ - min - / \xq{t + s) - yo{T)\ dr - S ] f{x). 

y s£[0,T] T Jo J 

Observe that 75 G C^([0, T] x R", K") is T-periodic with respect to time and 
so it can be extended from [0, T] to R by T-periodicity. By ([7]) we have now 



6 



O. Makarenkov, P. Nistri and D. Papini 



that 

Fs{9) = r \xo{r + e)-yo{T)\'dr- 
^' r (9) 



mm 



\xo{T + s)-yo{T)\'dT-TS. 



Under our assumptions the function 6 Fs(6) is T-periodic and contin- 
uously difFerentiable. Furthermore there exists a unique 6*0 G [0, T] such 
that Fo(0o) = 0. Since the function Fq reaches its minimum at 9o, we 
have that F^{9o) = 0. Then, there exists such that Fq{9) ^ for any 
9 G i9o- iyo,9o) U (6lo,6»o + i^a)- Thus F^{9) < for 9 e {9o ~ 1^0,^0) and 
Fq{9) > ioT 9 e {9o,9o + i^o)- Therefore for any S > sufficiently small 
there exists a unique 9s such that 

1) Fs{9s) = and F^iOs) > 0, 

2) 9s^9oasS ^ 0. 

Let > 0, by the previous considerations there exists (5^ > such that 



\xo{t + 0s) - yoMpdr - / \xo(t + 6*0) - t/oMpdr 

Jo 



<f (10) 



for any S G [0,(5^]. Then for a given S G [0,(5^], by applying the Malkin's 
results (Theorem 2 above) and by taking into account ([8]), we deduce the 
existence of a positive number es > such that, for any e G [0,6^], system 
(j4|) possesses a unique asymptotically stable T-periodic solution x^ such that 
Xfj.{t) G Ms^^xq), for any t G [0,T], satisfying 



\x^{t) ~ zjo(r)pdT - / \xo{t + 0s) - yo{T)\^dT 
Jo 



Finally, (fTU]) and pTjl conclude the proof. □ 

Remark 1 /n the case when equation has k solutions on the inter- 
val [0, T] we can reformulate Theorem 1 by simply replacing the sentence 
"unique asymptotically stable T-periodic solution" by "k asymptotically sta- 
ble T-periodic solutions" . Furthermore, we can interchange the role of the 
slave and master system in all of the previous arguments. 

3 A static feedback for the synchronization of 
trajectories 

In this section we consider two nonlinear nonautonomous systems 

i = 0(t,a;) (12) 

y = Ht,y), (13) 

where 0, V' : [0, +00) x R" — + R" satisfy the following conditions: 
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(3.1) t (j){t, x) and t — > x) are Lebesgue measurable functions for any 
X e R"; X —> 4){t, x) and x — > ^{t, x) are locally Lipschitz functions for 
almost all (a. a.) t > 0. 

(3.2) For any p > there exists an integrablc function 7p(t), t > 0, such that 

W,x)\<-fp{t) and |^(i,a;)| <7p(t) 
for a. a. t>{) and any x £ R" such that |a;| < p. 

(3.3) Any local solution x = x{t) and y — y{t) of the Cauchy problems 

( x = (t){t,x) I y = '^Kt,y) 

\ x{0) ^xoe R" \ y{0) = yo e R", 

can be extended to the interval [0, +oo). 

Let yo{t),t > 0, be the prescribed bounded solution of (fT3|) to which system 
([T2|) must be synchronized by means of a suitable feedback coupling term to 
be added to (|12p . The coupling term we propose here has the form of a static 
feedback given by 

Bsgn{x-yo{t)), (14) 

where B = diag(&i), hi < for any i = 1,2, and sgn(a;) = (sgn(xi), . . . , 
sgn(a;n)) for all x = (a;i)"=i £ R" . Therefore system (fT2|) takes the form 

X = (t){t,x) + Bsgn{x - yo{t)). (15) 

Since the right hand side of system (|15p is discontinuous in the state vari- 
able X we must adopt a suitable notion of solution for any Cauchy problem 
associated to (fT5|) . Here, following [9], as a solution of (fTSl) we intend an 
absolutely continuous function x{t), t > 0, such that 

x{t) e (f>{t, x{t)) + B K[sgn{x{t) - yo(i))] 

for a. a. t > 0, where 

K[sgn(a; - y)] = fl f| co (sgn(i3,(a; - y) \ A^)) , 

N C R" is an arbitrary set of Lebesgue measure zero, co{A) denotes the 
closure of the convex hull of the set A and Be (c) is the open ball in R" with 
radius e and center c. It is immediate to see that if to is a discontinuity point 
for some component of the vector sgn{x{t) — yo{t)), then the corresponding 
component of K[sgn(a;(io) — yo{to))] is the interval [—1, 1]. The other com- 
ponent of K[sgn(a;(io) — yo{to))] being -|-1 or —1. In the sequel we will refer 
to such solution as a Filippov solution to system p^ . 
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Then we associate to system ([T5|) the Lyapunov function V : R" :— 
{xgM:x>0} defined as follows 

n 
i=l 

where e = (e^)"^]^ is the tracking error given by Ci = Xi — yi, i — 1,2, n. 
We can prove the following result. 

Theorem 3 For any bounded set I C M" there exists a positive constant 
Ml > such that, if bi < —Mi for any i = l,2,...,n, then any Filippov 
solution x{t), t > 0, to system (llSp starting from a point xq ^ I converges in 

finite time to < ^— where fij — max{Af/ + bi, i — 1,2, ...,n} < 0, to 

2/0 (0 and x{t) = yo{t) for any t>tQ. 

Proof. Let x{t), i > 0, be a solution to (fT5l) with x{Q) — xq, whenever xq G /. 
Let e{t) = {ei{t))^^i , ei{t) = Xi{t) — yo,i{t), t > 0. From the chain rule for 
locally Lipschitz regular maps ([7], Theorem 2.3.9-(iii)), for a. a. i > 0, we 
have that 

lv{e{t)) = {C,m), 

for any ^ G dV{e{t)), where dV(e{t)) is the Clarke generalized gradient of 
the function V evaluated at e(t) ([7j, p. 27). In this case it can be easily 
seen that the vector ^ = is given by £,i — +1, (^^ = —1), for those 

indexes i G J+{t), (J-(t)), for which ei{t) > 0, {ei{t) < 0), while £_i G [-1, 1] 
for the indexes i G Jo{t) for which ei{t) = 0. Therefore, for a. a. < > such 
that e{t) ^ 0, we have that 

Jt^«t))= E e.(Osgne,(i) 

ieJ+{t)uJ-(t) 

= ^ [(j)^{t,x{t)) - ^/j^{t,yo{t)) +b,sgne,{t)]sgne^{t) 

ieJ+(t)uJ-(t) (16) 

= 53 [<j)z{t,x{t)) - ^p,{t,yo{,t))]sgne^{t) +b, 

ieJ+(t)uJ-(t) 
< max{M/ + fe;, z = 1, 2, n} = /i/, 

where Mi > \(j)i{t, x{t)) — ^i{t, yo{t))\ for a. a. t >0 and for any i = 1, 2, ..,n. 
Observe that, under our assumptions on the vector fields (j){t, x) and ^p{t, x), 
such a constant Mj > does exist. In fact, yo{t), f > 0, is a bounded 
trajectory in R" and it is not hard to see that any solution x{t), i > 0, to 
([T5|) starting from the bounded set / C K" is also bounded in M". By our 
assumption /i/ < 0, therefore by integrating (|16p between and t > we 
obtain 

V{eit))-V{e{0))<^IIt. 
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In conclusion, for t > > we have that V(e(t)) = and thus 

4 A dynamic feedback for the synchronization 
of trajectories 

In practice, since the switching is not instantaneous and the tracking error is 
not perfectly calculated, we have a serious drawback by using the previous 
approach, that is the so-called chattering phenomenon as shown in Example 
[2]of Section[3l To avoid this phenomenon we propose in the sequel a different 
coupling term to add to (|12p in order to solve our synchronization problem, 
namely a dynamic feedback coupling term defined by means of the singular 
perturbation theory. This method to eliminate the chattering and to make 
a controlled dynamical system insensitive with respect to external pertur- 
bations has been introduced in [4], where applications to specific tracking 
problems have been also presented. Specifically, given the reference trajec- 
tory yo{t), t > 0, which is a solution of system p^ . we introduce the function 
s : M+ X R" X M" R" as follows 

s{t,^o,x) = e^*(eo - 2/0 (0)) ~{x- yo{t)), 

where C is a given real symmetric matrix such that maxA(C) < —a < 0, 
here A(C) denotes the set of the eigenvalues of the matrix C. Then, under 
conditions (3.1) — (3.3), for e > small we consider the system 

{X = 0(<, x) — Bu 
ds ds (17) 

-^(^,^0,2^) + ^(i,Co,a:;)[0(^,a:) - Bu] := g(t,Co,a;,u). 

We are now in the position to prove the following result. 

Theorem 4 Assume, that B is a real negative defined n x n matrix. Then 
for any S > there exists e > such that 

lim \x,it) - yoit)\ < S, 

t — *oo 

where x^(t), i > 0, is the solution of |J7| ) such that ^^(O) = Re- 
proof. Let £ = and consider the corresponding system 
X =^ (f)(t, x) — Bu 

= Ce^'i^o - yo(0)) + yoit) - Ht, + Bu, 

Resolving the second algebraic equation with respect to u and substituting 
in the first equation we obtain 

xoit)~yait)^Ce^''i(o-yom, t>0 
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or equivalently 

2;oW-yoW = e^*(Co-yo(o)), t>o, 

with xo(t) satisfying xo{0) = ^o- By our assumption on the matrix C we 
have that hmt^oo(a;o(t) — yo{t)) = 0. Now we show that for given 6 > there 
exists e > such that 

\Xe{t) - Xo{t)\ <S 

for any t > 0. For this, we use the singular perturbation theory on unbounded 
time intervals. Specifically, in the sequel we verify that all the assumptions 
of the Theorem at p. 523 of (THI are satisfied for system (fT7|) in order to 
conclude that 

lim X!;{t) — xo{t) 

e—fO 

uniformly in [0, +cxd). First of all, for any {t,x) G x M" the equilibrium 
point 

u{ix) = -B-^ [Ce^*(eo - 2/0(0)) +V^(i,2/o(t)) - x)] 

of the equation edu/dr = g{t,^o,x,u) is asymptotically stable. In other 
words, the solution z = z{t), r > 0, of the Cauchy problem (the boundary 
layer) 

( z = g(i,£,Q,x,z) 
I z(0) = zo. 

converges asymptotically to u(t,x). In fact, consider the Lyapunov function 
M+ X K" X M" M+ given by 

Vii,x,z)^\^^^f^. 

We have, that 

^V{i,x,z{T)) = l^^V{i,x,z{T)),g{i,iQ,x,z{T))^ 

= {z(t) - u{i,x),g{i,^o,x,z{T)) - g(i,£,Q,x,u{i,x))) 
- {z{T)-u{i,x),B{ziT)-uii,x))) 
< -i/|z(t) - u(i,x)|2 

for some i/ > 0, since B is negative defined. Moreover, observe that the 
asymptotic stability is exponential and uniform with respect to (i, x) when 
X belongs to bounded sets. 

Furthermore, the origin a; = is an uniformly stable equilibrium point of 
(fT7|) when e = 0. To show this it is sufficient to observe that for e = we 
have 

x{t)-yo{t)^Ce'''{^o-yo{0)), t>0 
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and the change of variable e = x — yo makes the origin e = of this dynamics 
(exponentially) asymptotically stable, in fact 

(e,e) = (e,Ce) < -a|ep. 

At this point all the assumptions of the Theorem at p. 523 of [TB] are satisfied, 
thus 

lim Xs{t) — xo{t) 

E—>0 

uniformly in [0, +oo) and the conclusion of the theorem easily follows. □ 

Remark 2 The convergence theorem of f_?6] / employed in the proof of Theo- 
rem 4 also establishes that the absolutely continuous function u^, s > 0, given 
by (flT)) is such that 

lim Ue{t) = uo{t) 

£—►0 

uniformly on any interval [ti,+oo), ti > 0, where 

uo{t) = [Ce^*(Co - yo(0)) + 2/0 W) - <^(i, -^oW)] , t>0, 

is the so-called equivalent control of the theory of variable structure systems 
(see \23 ] and 124-1)- In other words, it is the ideal control which realizes 
exactly the condition s(t,xo(t)) ~ 0. This convergence property together with 
the absolute continuity of the functions and uq prevents the chattering 
phenomenon. 

Finally, we would like to point out that the proposed dynamical feedback 
coupling term —Bu, designed by means of (fT7)l . makes the dynamical system 
X = (j>{t, x) insensitive with respect to (possibly unknown) perturbations. 
In fact, consider any bounded perturbation p : R+ x M" — > M" of system 
X — (t>{t,x) satisfying (3.1) — (3.2). Thus system P7|) takes the form 

{X = 4'{t,x) -\- p{t,x) — Bu 
ds , ^ , ds , ^ , , , ^ , (18) 

^ -^(i^Sca;) + — (t,^o,a;)[<P(^,a;) +p(i,a;) - -BuJ. 

It is easy to verify that for this system we can repeat the arguments used 
in the proof of Theorem 3. In fact, for e = 0, the dynamics of the error 
e = X — yo for system ^TE\\ is the same of that of system P7)) , namely 

i(i)-yo(0=C^e^*(eo-|/o(0)), i > 0, 



Moreover, u{i,x) is still an asymptotically stable equilibrium point of the 
second equation of system 
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5 Examples 



In the following examples we have chosen to show the behavior of only the 
first component of the solutions of the slave and master systems, the other 
components having a quite similar behavior. 

Example 1 At first we illustrate the result of Section [H on the phase syn- 
chronization in M^. We take the same system as slave and master systems, 
namely 

1 , 9 



Xi 



Xi 



X2 



20 



Xi 



7 
10 



(19) 



In 



which has an asymptotically stable limit cycle with period T ~ 9.83 
particular, let yo be the T -periodic solution of (jl9p which starts at t — 
from the point (—0.7481, 1.5164) of the limit cycle. In this way, since 

l-T 

system ^ becomes 



1 



mm 7^ / |?/o(t + s) -?/o(T)pdr = 0, 

se[o,T] T Jq 



Xi 



X2 



2 xi 



1 3 



X2 



-1(1 

20,' ^ 



llxi 



^ e{\x ~ y,{t)\^ ~ 5) 
\x-yo{t)\^~S) 



(20) 



and we consider the solution x of (j20p which starts from x{0) = (5,-5), 
with the choice e = 0.01 and 6 = 0.05. In Figure Q] the solid line is the 
graph of the first component x\ of x in the intervals [T,2T] (left picture) 
and [SOT, 5 IT] (right picture), while the dashed line is the graph of the first 
component of the T -periodic solution yo- Due to the asymptotical stability 
of the limit cycle of (|19p . we see that xi already has the shape of the first 
component of yo after only one period, but the phases of these trajectories 
are still considerably different. The picture on the right shows that the phase 
difference between the first components of x and yo is significantly reduced 
after 50 periods thanks to the coupling term added in (j20p . 

Example 2 We exploit the result of Section\^by considering the static feed- 
back synchronization of the following neural network introduced in \29i/ 



(21) 



where f{s) = (|s + 1| — |s — l|)/2. In 129^ it is shown that system (f2T|) has a 
chaotic attractor. As master system we use another neural network 
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+ I{t), (22) 
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where I{t) is the periodic input such that (j22p has the 2tt -periodic solution 
yo(t) = (cos i, sini, — cosf). We take the gains 61 = 62 = ^3 = —3-5 in 
the coupling term (|14p and consider the solution x of (|15p such that x{Q) = 
(— 1, 1, 1); then in Figure\^ the graph of the first component xi of x is plotted. 
In particular, the picture on the left shows the convergence in finite time to 
the first component of the reference solution z/q; on the other hand, the picture 
on the right is a zoom of the left one around the hitting zone and is a clear 
evidence of the chattering phenomenon. 

Example 3 To illustrate the dynamic feedback synchronization considered in 
Section and to compare the issues with those of the previous example, we 
consider again the two neural networks (PT|) and of Example with the 
same forcing term I{t). We take C — B = diag(— 1, —1, —1) and e = 0.001 
and consider the solution x of (fT7|) such that x{0) — (—1, 1, 1). In Figure\^the 
graphs of the first components xi of x (solid line) and cost ofyo (dotted line) 
are plotted. In particular, the picture on the left shows how xi approaches 
cos t, while the picture on the right is a zoom of the previous one around 
t = 2iT and shows that xi{t) remains in a neighborhood of cost without any 
chattering. 

Remark 3 Observe that in Examples\^and\^we have taken the initial con- 
ditions for the master system on the reference trajectory yg. On the other 
hand, in | _/<$] / it is shown that yo is globally exponentially stable, therefore 
these initial conditions could be also chosen as far as we like from yo- 




Fi gure ll Time behavior of xi[t) (solid line) in Example [Tl in the intervals [T, 2T] (left 
picture) and [SOT, SIT] (right picture). The dashed line is the graph of the first component 
of the T-periodic solution yo{t) of the master system. 
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Figure 3: Time behavior of xi(t) (solid line) of Example |3] The dotted line is the graph 
of cos t, that is the first component of the periodic solution of the master system. 
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